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Abstract 

We consider the large N limit of three-dimensional U{N)j t Chern-Simons 
theory coupled to a Dirac fermion in the fundamental representation. In this 
limit, we compute several correlators to all orders in the 't Hooft coupling 
A = N/k. It was suggested recently that this theory is dual to the Legendre- 
transformed theory of scalar fields coupled to Chern-Simons gauge interac- 
tions. Our results show that this duality holds for any value of the 't Hooft 
coupling, at least at the level of the planar 3-point functions. In addition, 
we determine the sign in the duality transformation of the Chern-Simons 
level, as well as the relation between the "triple-trace" deformation which 
exists in the bosonic Chern-Simons theory and in the Legendre-transformcd 
fermionic theory. 
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1 Introduction 

Let us consider three-dimensional free field theories with N massless complex 
bosons cf) 1 or Dirac fermions These conformal field theories are known as the 
"free U(N) vector models" , and we will restrict ourselves to their U(N) singlet sec- 
tor. The closely-related "critical U(N) vector models" are interacting conformal 
field theories that can be defined (in the large N limit) by performing a Legendre 
transform of the free theories with respect to the scalar primary operators or 
i[) l ipi. Because we gauged the U(N) symmetry, these four theories can be naturally 
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coupled to a gauge field with Chern-Simons interactions at level k. The singlet 
sector of the vector models is recovered by taking k to infinity. 

In the large N limit with fixed 't Hooft coupling A = N/k, these theories have 
remarkably simple properties. First, both the bosonic and fermionic theories are 
conformal 1 . Second, the full spectrum of gauge- invariant local operators in these 
theories is independent of A at large iV [1, 2]. The structure of the planar 3- 
point functions of these operators is also highly constrained [3]. In particular, [ ] 
determined all the planar 3-point functions of primary operators in terms of two 
parameters 2 N and A, which are some functions of N and A. 

The results of [3] also imply that the planar 3-point functions of the free 
fermionic vector model turn into those of the critical bosonic model, when one 
changes the coupling A from zero to infinity. A similar relation holds between 
the free bosonic and critical fermionic vector models. This three dimensional 
"bosonization" was first conjectured in [ ] by considering the gravity duals of 
these theories. In particular, it implies that (at least at the level of the planar 
3-point functions) the U(Nf)k f fermionic Chern-Simons vector model is equivalent 
to U(Nb)k b critical-bosonic Chern-Simons vector model, for some integers Nf,kf 
and N b ,kb- 

In [4], the parameters N and A of the bosonic Chern-Simons theories were 
determined exactly as functions of iV and A. Using this result it was found that, 
in the large N limit, the U(N)k critical-bosonic Chern-Simons theory is equivalent 
to the U(\k\ — iV)±fc fermionic theory. Note that the pure Chern-Simons theories 
on both sides of this relation are equivalent by level-rank duality [5, (3, 7], which 
is valid for any iV and k. The duality takes a more familiar form when written 
in terms of a shifted level /cym, defined by k = &ym + sign(fcYM)^; it is given by 
U(N)k YM U(\kyu\)N (up to the level's sign), where the subscript denotes the 
value of the shifted level. This is an important consistency check on the proposed 
duality in the large N limit; in some computations, such as correlation functions 
of Wilson lines, the leading contributions are of order N 2 and are completely 

1 At finite N the fermionic theory is still conformal. The bosonic theory has a possible defor- 
mation (cj) l (j)i) 3 that is marginal in the planar limit, but not at finite N. At weak coupling one 
may nevertheless reach a fixed point by tuning the coupling of this deformation [ ]. 

2 In the bosonic and critical-fermionic Chern-Simons theories there is a third parameter, cor- 
responding to the deformation mentioned above, that only affects the 3-point function of the 
primary scalar operator. 
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determined by the pure Chern-Simons contribution. 

The duality relation above was found by taking the strong coupling limit of the 
exact expressions for N and A in the critical bosonic theory, and comparing the 
result to the same parameters in the weakly coupled fermionic theory (which were 
determined perturbatively in [2]). The relative sign between the Chern-Simons 
levels in this relation was not determined. In this work we repeat the computations 
of [4] for the case of the fermionic Chern-Simons vector model. We find that 
the above duality holds for any value of the coupling (and not just in the weak 
coupling limit of the fermionic theory), and we determine the relative signs between 
the levels in the duality relation. Moreover, the bosonic (and critical-fermionic) 
theories contain a "triple-trace" deformation, which in the bosonic theory can be 
written as Ag(0 J 0j) 3 . This deformation is exactly marginal at large N [1], and we 
determine how it maps between the two theories. 

Let us summarize our results. We find that the planar 3-point functions of the 
U(N)k Chern-Simons theory coupled to a fundamental scalar field in the critical 
fixed point, are equal to the same correlators in the U(\k\ — iV)_fc Chern-Simons 
theory coupled to a fundamental fermion. This provides further evidence that 
these two theories are equivalent 3 . The same duality map exists between the 3- 
point functions of the bosonic and critical fermionic Chern-Simons theories. In 
addition, the "triple-trace" deformations of the latter two theories are related by 

A* = 8^(1 - |A/|) a (3 - SttA/AJ) . (1) 

The paper is organized as follows. In section 2 we introduce our conventions and 
methods. In sections 3 and 4 we present the detailed computations of correlation 
functions in the fermionic theory. The reader interested in the final results can 
skip directly to section 5 in which we use our results to determine the duality map 
between the fermion and boson theories. In section 6 we present the results for the 
critical fermionic theory. We conclude with a discussion in section 7. Appendix A 
includes a discussion of reflection positivity in theories with a non-real action. 

3 At finite N, the conjectured duality must involve a half-integer shift of the Chern-Simons 
level of the fermionic theory due to the parity anomaly [4] . 
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2 Fermionic Chern-Simons Vector Model 



Consider the theory of a Dirac fermion ip in the fundamental representation of 
U(N), coupled to gauge bosons with Chern-Simons interactions at level k in 
three Euclidean dimensions. The action is 4 

S = -tj d ' Xe " UP ( A ^p + \f abc AlA b v A*\ + J d 3 x^V^ . (2) 

We work in the 't Hooft large N limit, keeping A = ^ fixed. The normalization 
of the Chern-Simons action (2) is such that the theory is gauge invariant if k G Z, 
up to a ±| shift due to the parity anomaly [8, 9, 10]; this shift will not matter 
to us since we work in the large k limit. We use light-cone gauge, A_ = 0, where 
the light-cone coordinates are defined by x ± = x T = (x 1 ±ix 2 )j^2. As was 
first noticed in [2], this gauge choice leads to great simplifications in perturbative 
calculations, since the self-interaction of the gauge bosons vanishes. 

In this gauge the gluon propagator, which is exact in the planar limit, evaluates 

to 

(A;(-p)Al(q)) = G v ,{p)5 ab ■ (2n) 3 5 3 (q - p) , 

G +3 (P) = -G 3+ (P) = 4 ^^, (3) 

k p + 

and the other components of G vanish. The exact planar fermion propagator in 
light-cone gauge was computed in [2] by summing over rainbow diagrams, and was 
found to be 

^iipWi-q)) = 5jS(p) ■ (2n) 3 5 3 (p - q) , (4) 
S{P) = 2 , (5) 

where p 2 s = 2p + p~ = p\ + p\. 

As in [4] we use dimensional regularization in the x 3 direction and a cutoff A 
on the momenta in the radial direction of the x 1 — x 2 plane. This regularization, as 
well as the light-cone gauge-fixing condition, preserve an SO (2) rotation symmetry 

4 The covariant derivative is 2? M = + A®T a . The generators T a are anti-hermitian, in the 
fundamental representation of U(N), and normalized such that ti(T a T b ) = — \$ ah ■ The Dirac 
matrices are the ordinary Pauli matrices 7^ =a», fj,= 1,2,3, and ip = ip^ . 
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in the x l — x 2 plane. While this cutoff prescription breaks Lorentz invariance 
(as does our choice of gauge), conformal invariance and gauge invariance, these 
symmetries are restored by demanding that the theory is conformally-invariant in 
the continuum limit. 

The spectrum of primary operators of this theory contains a single primary 
operator for each integer spin s > 1, with dimension A = s + 1 + 0(1/N), and 
a scalar of dimension A = 2 + 0(1/N) [2]. All other primaries are obtained 
by products of these "single-trace" operators, and at large iV their dimension 
is just the sum of the dimensions of the individual single-trace operators in the 
product. In this paper we will compute correlation functions of the spin and 
spin 1 operators 

J (0) =^V,, (6) 
J« = i^fc . (7) 

As we explain in appendix A, the factor of i is included in to make its 2-point 
function positive-definite (as it is in the critical bosonic theory defined in [4]). 

3 Exact Vertices 

In this section we compute the exact vertices of J(°\ and J_ , which are 
defined by the correlators 

(J {0 \-q)Mk)fi(-p)) = V(q,p)5j ■ (27r) 3 5^\q +p-k), (8) 
(j£\-q)MW(-p)) = V±(q,p)6i ■ (2n) 3 5^(q + p - k) . (9) 

Here we suppress the spinor indices on fermion fields and on V and V±, which 
are matrices in spinor space. Moreover, in the above correlators we amputate the 
external fermion propagators. 

The computation is done by solving the large iV Schwinger-Dyson equations 
for V and V± in light-cone gauge. The equations are shown in Figure 1. We set 
q* 1 = 0, which leads to simplifications in the calculations. These exact vertices will 
be used in section 4 as building blocks for the gauge invariant correlation functions 
of J(°) and JW. 
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— (g) = — - X + 

p 

Figure 1: Schematic form of the Schwinger-Dyson equations for the vertices. The 
circled cross denotes an insertion of an exact or vertex. The un-circled 
cross denotes an insertion of the free vertex. A filled circle denotes the exact 
fermion propagator. 

3.1 J(°) Vertex 

Let us first consider the scalar operator vertex (8). The bootstrap equation in 
Figure 1, with the circled cross denoting V(q,p), is given by 

— ^S(k + q)V(q, k)S(khM-^—^ . (10) 

The vertex V(q,p) is a 2-by-2 matrix in spinor space. It can be expanded in a 
complete set of matrices consisting of 7 M and the 2x2 identity matrix / as 

V = v^ + vjl. (11) 

Note that a general 2-by-2 matrix A = a^^+ajl satisfies 7t 3 ly47' + ] = 2 (a/7 + — a_7). 
It then follows from (10) that t>_ = v% — 0. Plugging (11) and the fermion prop- 
agator (5) into (10), performing the Dirac algebra, and separating the / and 7 + 
components, we obtain 

f d 3 k Vl (q, k) {2i\k s - q 3 ) k+ + 2v + (q, fc)(fc+) 2 
Vl (q,p) = 1 - 4«A J ^ k 2 (k + q) 2 (k-p) + ' (12) 

d 3 k Vl (q, k) [(1 - 2\ 2 )k 2 + k 3 (k 3 + ? 3 )] + v+(q, k) (q 3 + 2i\k s ) k + 



(2tt) 3 k 2 (k + q) 2 (k -p)+ 

(13) 



From these equations and our choice of external momenta we can see that V(q,p) 
can depend only on p and q 3 . Furthermore, from dimensional analysis and rota- 
tional invariance in the x 1 — x 2 plane, V(q,p) can be written as 



V(q,p) = v+(q,p)-f + +Vi(q,p)I = 2 ^ ± g ( r^r ) 7 + + / \ri\ 1 

Ps \m\J \m\J 



(14) 
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where the factors are for later convenience. Our remaining task is to compute / 
and g. 

Let us plug the ansatz (14) in (12), (13). We perform the integral over fc 3 , and 
then over the polar angle in the k 1 — k 2 plane 5 . We are left with the integral 
equations 

2p.\ _ r\ t (»M > -<fiO/(i»)+2Mi5D 



/ ^ = l-2iA/ dk s -^- 2 (15) 



2X P+g /2pA = _2^ dk 2k * (1 - * 2 ) / (&) + (& + 2<A*.) (g. 



\M/ P + JO " ^ k s+Q3 

(16) 

It will be convenient to work with the dimensionless, parity-invariant 6 variables 
y = f 2 !, x = pH, A' = rn-, and A = A ■ sign(g 3 ). We obtain 

/(y) = 1 - zA fdx ^9(x)- (l- iXx)f(x) 



ry 

yg{y) = -* / 



1 + x i 

y , x 2 (l-A 2 )/(x) + A(l + zAx)x^(x) 



1 + x 2 

These equations can be solved, for instance by using Mathematica, and we find 

^ _|_ g— 2iA arctan(j/) 
fW' - i _|_ e -2iAarctan(A') ' 

_ i _ f Ay - (1 + iAy)e~ 2darctan ^ 

Ay ]_ _|_ g— 2iA arctan(A') 

3.2 Vertices 

The computation of the J± vertices (9) is similar to that of the vertex. The 
Schwinger-Dyson equation for = iip t '-f tJ "i(j i is 

/d 3 k 1 
— - 3 ^S(k + q)V»(q, fc)5(*) 7 l+] ■ (21) 



5 The polar coordinates k s ,8 are defined by fc ± = k s e l6 /y/2. 

6 Our theory is invariant under a parity transformation combined with A — » —A. 
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As before, we choose q* 1 = 0, and the vertices are independent of p%. Based on 
similar considerations as in the case, we can write V + as 



V + (q,p) = 9 {+) (yh + + 2 -^f {+) (y) 

Q 



(22) 



Plugging this in the bootstrap integral (21), and carrying out the k 3 and polar 
integrals as before, we find 



1 + X 2 



' x\l - \ 2 )f (+ \x) + (1 + i\x)g ( -+\x) 



dx 



1 + x 2 



(23) 
(24) 



The solution to these equations is 



f (+ \y) 



l-e 



2iA(arctan(A') — arctan(j/)) 



gW(y) = - l-i\y + (l + ^?/)e 2d{arctan(A ' ) - arctanfo)) 
Similarly, V~ can be written as 

V-(g,p) = ij- + 2 ^g^(yh + + 2 -^f^(y) . 

Pj 93 
The equations that result from (21) are 



(25) 
(26) 

(27) 



y 2 f(-\y)=i\ / dx 



y x 2 n ( 



'■g(-\x) - (1 - i\x)x 2 f(-\x) - z(l + \ 2 )x 2 + 2Xx 
1 + x 2 



(2f 



., ,_, f y (l-\ 2 )x i f(-\x) + (l + i\x)x 2 g(-\x)-i(l-\ 2 )x 2 (l-i\x) 

y 9 yy) — ^ / ^ 



1 + x 2 



and their solution is 

f { -\y) 
g { -\y) 



i 



1 - 2?Ay - e - 2darctan ^) 
(1 + zAy)e- 2darctan{y) - (1 - zAy) 



(29) 

(30) 
(31) 
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4 Correlation Functions 

Using the results of the previous section, in this section we compute the exact 
planar 2-point and 3-point correlators of and J^ 1 ' with the external momenta 
all in the x 3 direction. The 2-point functions can be computed by the diagram 
in Figure 2. Note that only a single insertion of an exact vertex is required in 
the 2-point function to account for all the Feynman diagrams without any double- 
counting. The 3-point functions are obtained by computing the diagrams shown 
in Figure 3. 




Figure 2: Diagrams contributing to (j( s >j( s >). 




Figure 3: Diagrams contributing to ( and ( J W ). 
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4.1 2-point Functions 

4.1.1 (J(°)J(°)> 

Using the exact vertex in (14), (19), (20) and the fermion propagator (5) the 
2-point function evaluates to 

(j(°)(- g )J<°)) = -N J^- 3 TT[S(q + p)V( q ,p)S(p)] 
Ntm(^). N A 

where the overall minus is from the fermion loop, and in the last step we took 
A — > oo, keeping divergent terms. 

Let us denote the background source of J(°> as if. The divergence in (32) can be 
subtracted with a mass counterterm f> 2 , whose coefficient is uniquely determined by 
demanding that the continuum theory be conformally-invariant. The renormalized 
2-point function is therefore 

(■7 (0) (-g)^° ) ) = - JVt Xp ) |g|- (33) 

This result is also correct for general q. In appendix A we verify that it is consistent 
with reflection positivity. 

4.1.2 (J«J«) 

Similarly, using the exact vertex (22), (25), (26) and the fermion propagator 
(5), we obtain 7 

= -Ni /tS^Tt [S(p + q)V + (q,p)S(p) 1 -] 

Niq(e m ' x -1) N k 

-»■ — + —A . 34 

16vrA 4tt v ; 

The divergence can be removed by adding a mass term to the background gauge 
field that couples to The A-odd part of (34) is 



8 vrA 

7 One could also use the exact vertex in this computation and get the same result 
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This parity-odd 8 piece can come from a conformally-invariant contact term (J fl (q)J u ) ~ 
e Atjyp g p , which corresponds to the appearance of a Chern- Simons term |^ J A A dA 
in the generating functional F[A, ...]. It was argued in [11] that the fractional 
part of this term (in units of k) is a physical observable. We will comment on this 
issue in section 7. 

Removing the divergence and the contact term, we are left with the parity-even 
result 

m _ q)J?) = _Nj^ M (36) 

This 2-point function is uniquely determined by conformal invariance and current 
conservation up to an overall coefficient, and we may therefore write 

( j(D ( _ ff) jJD) = . q ^~^ q2 . (37) 

In appendix A we verify that this result is consistent with reflection positivity. 

4.2 3-point Functions 
4.2.1 (J(°)JWJ(°)> 

Evaluating the diagrams shown in Figure 3 (a) using the exact vertex, we 
obtain 

(J i0) (-q)J {0) (-qV i0) (-Q")) = ~N J ' ^Tr[s(p + q)V(q,p)S(p)V(q" ,p)S(p - q")V(q' ,p) 

+ S(p + q)V(q,p)S(p)V(q',p)S(p - q')V(q",p) 

N tan 2 (^) t 
= — 2 J . (38) 

Z7TA 

This is a pure contact term. For A = this correlator vanishes since 

J(0) 

IS 

parity-odd. For non-zero A the vanishing of this correlator was predicted in [ 
by computations in Vasiliev's gravity theory and by a perturbative calculation in 
field theory. This result is also consistent with the general analysis of Maldacena 
and Zhiboedov [31. 



J Here we are referring to pure parity transformation, without the additional A — > —A. 
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4.2.2 (jWjfjW) 
Finally, using the and J±^ vertices in Figure 3 (b), we obtain 



(J {0 \-q)J?(-qV { -\-q")) = ~N 



d 3 p 



Tr 



S(p + q)V(q,p)S(p)V-(q',p)S(p-q')V + (q",p) 
+ S(p + q)V(q,p)S(p)V + (q",p)S(p - q")V~(q' } p) 



N 1 q" 2 [1 - e^ +X ")] + q' 2 [e wi ' x - e"^'] 



47rA 1 + e 7 ™ 



q'q" 



(39) 



This correlator has one parity-even and one parity-odd structure, which are conformally- 
invariant [12]. The A-even part (which corresponds to the parity-odd structure) 
is 



(J (0) hq)J { +\-qV ( -\-q"))x 



and the A-odd part is 



(J i0 \-q)J?\-q')J { -\-q"))x-oM 



iN sin(7rA) 
8^A 



111 - K 
q" q 



//i /ll 
q ' 



q' q" 



(40) 



N sin 2 



4ttA 



, , q t q , q 

q'q" qq" qq' 



(41) 



The last term in the brackets is a contact term that we discard. 



5 Analysis of the Results 

Let us summarize our results, omitting contact terms. First, the 2-point functions. 



Fermion 2-point functions: 

(J(O)J(O)) 



7rA 

2iVsin(7rA) 
7rA 



(J(O)J(O)) 



fcr. 



fcr. 



(42) 
(43) 
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Here, (-)f er . are equal to one half times the same correlators in the theory of a 
single free Dirac fermion. They are obtained by setting A = and N = 1/2 in the 
results of section 4.1. Next, we write the normalized 3-point functions, defined by 



Fermion 3-point functions: 



(44) 



(j(°)j(Vj(°)y 

(jmjWjWy 



(45) 



ttA 



4iVtan(^) 
+ sign(A) 



1/2 



(j(o)j(i)j(D) 

1/2 



norm, 
fer. 



7T 



Atan(^) 



AN 



(j(o)j(i)j(D) 



norm, 
odd 



(46) 



Here, (■}fe° ri11 ' * s again half the normalized correlator in the theory of a free Dirac 
fermion, and OodT 1 ' ^ s an independent structure that only appears in interacting 
theories. In fact jW j(O) 

odd, for s, s' > 0, is equal to the single conformally- 
invariant structure which appears in the same correlator in the critical bosonic 
vector model. We therefore define this structure to be half the normalized corre- 
lator in the theory of a single critical complex scalar. This last correlator can be 
obtained by taking A = and iV = 1/2 in the result for critical bosonic Chern- 
Simons vector model, computed in [4]. 

As in the bosonic theory [ ], we can use these results to determine the pa- 
rameters N and A of Maldacena and Zhiboedov [3] 9 . Since J^ 1 ' is canonically 
normalized 10 , we can read off N from the J*- 1 -* 2-point function. Then, A can be 
computed from (j(o) j(i) j(i))norm._ We find 



N = 2N 



sin(7rA) 
ttA 



A 



/ttA 
tan — 
V 2 



(47) 



9 In [3] only correlation functions of even-spin operators were computed, and we assume (as 
in [4]) that these results extend naturally to a theory with operators of any integer spin. 

10 The current J 1 - 1 ' is canonically normalized in the sense that [Q,i/j] = ifj, where Q is the U(l) 
charge defined by the current. 
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5.1 Comparison with Critical-Boson Theory 

We now determine the duality relation between the fermionic and critical bosonic 
theories. In the fermionic theory, let us denote A by A/ and iV by Nf. In [4], 
the variables A and N were computed in terms of the parameters A& and Nf, of 
the critical bosonic theory. They are given by N = 2iVbsin(7rAfe)/7rA{, and A = 
— cot(7rAb/2). Note that this A has a minus sign compared with the one written 
down in [ i] . It comes from matching the conformal structures of the critical bosonic 
theory to the ones we are using in the fermionic theory 11 . 

We equate the parameters of the critical bosonic theory with (47), 

„ T sin(7rAf) , T sin^At) . , 

N = 2N f — = 2N b — , (48) 

TlXf TTAf, 

X = ta.££) =-«*££). (49) 

Notice that the 2-point functions of in the two theories change sign at A/, A& = 
±1. Therefore, the theories are unitary only when |A/|, |A&| < 1. Now, equation 
(49) implies that A/, A& have opposite signs, and that 

- »(^)=-»(^). <*>) 

The solution is |A/| + \X b \ = 1, and we also see that sin(7rA/) = — sin(7rAb). Using 
A = N/k, equation (48) then gives 

k f = -k b . (51) 

Finally, we have 

N f = \k b \-N b . (52) 



n To see this, one can use the results of [3] (and the values of N and A written above) to compute 
the even structure in ( jW This is donc by first computing ( j(°> 
where — \ is a parity-odd scalar operator. This is necessary because in [ ] it is assumed 
that the dimension 2 scalar operator has odd parity. The result can then be verified to agree 
with the direct calculation of [4], eq. (70). 
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6 The Critical Model 



In this section we consider the critical fixed point of the fermion theory, the 3D 
Gross-Neveu model with Chern-Simons interactions. We will define this theory 
in the planar limit as the Legendre transform of the theory (2) with respect to 
the scalar operator ipip. 12 We implement this by introducing a coupling aipip, and 
adding a path integral over a. The spectrum of the critical theory at large N is 
identical to that of the regular theory, except that the scalar operator (which is 
now a) has dimension 1 + 0(1/N). 

For the sake of clarity, let us briefly review how to compute correlators in the 
critical theory Consider first the generating functional of connected correlators in 
the regular theory (let us call it S e s.), that is defined by 

e -s cB . (a,A) = JvipVip VA e - s -f d3x °i>i>-f d3 * -M . (53) 

Here, A collectively denotes sources of currents with s > 0. In previous 
sections we computed S e ft. at large N, up to cubic order in the sources. 

To define the critical theory we make a dynamical by including a path integral 
over it, so that S e g_ becomes the effective action for a. In the large N limit - 
where we include only planar contributions on the right-hand side of (53) — the 
effective action is proportional to N, and it is therefore the c-lPI effective action 
of the theory 13 . For example, the 2-point function of a is given by the inverse of 
the quadratic o term in S e s., namely 

G(q) = (a(-q)a) CTit . = (-{jMj^U^.)- 1 = ^^^ l . (54) 

The other 2-point functions, (j( s 'j( s >) with s > 0, are unchanged in the critical 
theory. 

Before proceeding, we note that the critical theory admits a cr 3 deformation 
that is exactly marginal at large N. It is analogous to the "triple-trace" (0*0;) 3 

12 One way to think of the critical theory is as the UV fixed point at the end of an RG flow that 
starts with the regular theory, deformed by the operator (tpip) 2 . However, since this deformation 
is irrelevant, it is difficult to make this description into a rigorous definition. 

13 Beyond the large N limit one has to include 1 /N corrections to S e g, from non-planar diagrams 
in the regular theory. In addition, the 1PI action of a now receives contributions from a loops, 
so the critical and regular theories are no longer related by a Legendre transform. 



16 



deformation of the bosonic theory. The effect of turning on this deformation is to 
change the effective action by 5S e {{. = N J d 3 x |r°" 3 - At large N, this interaction 
only affects the a 3-point function. 

To compute general connected correlators that involve a, one may take the 
Legendre transform of S c s, with respect to a; the result is the generating functional 
of the critical theory. Up to 3-point functions, we may equivalently compute a 
connected correlator by reading off the c-lPI correlator from S e g,, and multiplying 
each a leg by its 2-point function (54). 

We are now ready to compute the parameters iV q b and A q .b. of the critical 
model, which is a "quasi-bosonic" theory in the language of [3]. Since the 2-point 
function (jWjW) i s unchanged, the result (47) for N still holds, and we have 

sin(7rA) . , 

iV q . b . = 2iV^-^ . (55) 

TV A 

To extract A q .b. we follow the conventions of [3] and consider a scalar operator 
with even parity by defining a = NX q ±,a. The factor of N is added so that a 
correlators will be of order N in the planar limit. Using (54) and (41), we find 
that 



Irr /( 1 )\ norm - — s ^S n (Vb.A) 

\U J J /crit. fer., even r> 



vrAtan (^) 



N 



1/2 

/ 7(0) 7(1) 7(1)\ 

free bos. 



(J (0) J (1) J (1) )™' U 



(56) 



On the left-hand side we have the parity-even structure of the critical fermion 
model. The correlator on the right-hand side is one-half the normalized 3-point 
function in the theory of a free complex scalar. Using the results of [3], one can 
see that the right-hand side of this equation must have a positive coefficient, and 
therefore sign(A q .b.) = — sign(A). Completing the calculation as in section 5, we 
find that 

Aq. b . = -cot^J . (57) 

These results mirror those of the critical bosonic model [4]. In terms of k and N, 
the duality map is exactly the same as the one relating the regular fermion and 
the critical boson, appearing in (51) and (52). 
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Finally, let us determine the duality map for the triple-trace coupling 14 . To 
this end we compare the 3-point function of a with that of in the regular 
bosonic theory. For the 3-point function of a, using (38), (54) and (57) we obtain 

(a(-q)a(-q')a(-q")) = N%M-q)a(-q')a(-q")) 

= N 3 \l h G(q)G(q')G(q") \-NX f 6 - (j(°)(-q)j(V(-q>)j(°)(-q»)} 

1 



32iV7r 2 A 2 cot 



VTA 

T 



1 + 27tAA 6 cot" 



VTA 

T 



9 9 9 



(5? 



Normalizing this correlator by the 2-point function (54), equating to the corre- 
sponding bosonic correlator (j(o) j(o) j(o)^norm. that 

was computed in [4], and ap- 
plying the duality map for k and N, we obtain the relation between the triple-trace 
couplings in the two theories, 



A 



8vr 2 (l - |A/|) 2 (3-SttA/AJ) . 



(59) 



7 Conclusions 

In this work we obtained all the planar 3-point functions in the fermionic Chern- 
Simons vector models, by computing several correlators exactly and using the 
results of [3]. Using this result we completed the duality map of these theories to 
the bosonic Chern-Simons vector models, and verified that it works for all values 
of the 't Hooft coupling and of the triple-trace coupling. 

Since we worked in momentum space, the results of our computation include 
the contact terms that appear in the 2- and 3-point functions. Contact terms 
are usually not universal, but in [J 1] it was noted that certain contact terms, 
specifically those in (J«J«) and (J^J®), are in fact physical observables (up 
to integer shifts). One can therefore use these terms to further test the duality. It 
would also be interesting to understand if contact terms in other 2-point functions 
are physical, and whether they are constrained by high-spin symmetries. Moreover, 
the coefficients of these terms map under the AdS /CFT correspondence to the bulk 
0-angle and its higher-spin analogs in Vasiliev's theory. Since these terms do not 

14 We are grateful to S. Giombi for reminding us of this point. 
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affect the bulk equations of motion, they were not known previously, and it will 
be interesting to analyze them. 

There are several issues one has to take into account when trying to match 
these contact terms under the duality between vector models 15 . First, the gauge 
and topological £7(1) currents can mix and one has to determine the combinations 
which map to each other under the duality. While this effect cannot influence 
correlators at separated points, it will in general affect the contact terms. Note that 
the contact term in (J^J^) does not suffer from such an ambiguity. Moreover, 
since our gauge and regulator break Lorentz invariance the correlators may contain 
Lorentz- violating contact terms in addition to the ones written down in [11]. 

It is currently not known how the results of [3] extend to supersymmetric 
theories. It will be interesting to understand this, both using the current algebra 
techniques of [•>], and by directly computing the correlators as we did here. For 
instance, it will be interesting to see how this duality plays out in theories with M = 
1 supersymmetry, as there are no similar known dualities in this case. We note that 
the form of the "bosonization" duality obtained in [ 1] is quite similar to Giveon- 
Kutasov duality [13] of M = 2 supersymmetric theories, or one of its "crural" 
versions [1 J]. It would be interesting to understand if the non-supersymmetric 
duality of this paper follows from the supersymmetric one. Analyzing M = 6 
theories in these methods would also be interesting, in particular given the relation 
between Vasiliev's theory and string theory proposed in [15]. 

So far the tests of this duality mostly follow from the large N high-spin symme- 
tries of Chern-Simons vector models 16 . It will therefore be interesting to find more 
evidence for the bosonization duality which does not follow only from symmetry 
considerations. One direction is to attempt to perform large N computations of 
partition functions on various manifolds, or to consider correlators of Wilson loops, 
as these are not entirely determined from symmetries as far as we know. 

One such possibility is to compute the thermal free energy at large N. This 
was done in [2, 16], but the results were found to be inconsistent with existing 
(conjectured) dualities. One possible reason for the discrepancy, already mentioned 

15 We thank Zohar Komargodski and Nathan Seiberg for discussions on this issue. 

16 The fact that the duality transformation at large N is the same as level-rank duality, and 
the matching of mass deformations done in [4], do not follow from the symmetries as far as we 
know. 
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in [16, 4], is that the holonomy of the gauge field along the thermal cycle might be 
non-trivial (as assumed in these computations). In [17] it is shown that, indeed, 
when one carefully accounts for the holonomy the dualities hold. 

It will be very interesting to find evidence for the duality beyond the large N 
limit. At this point, it is not yet clear what are the high-spin symmetry constraints 
on non-planar correlators, or on planar n-point functions with n > 3. One may 
try to understand whether such constraints exist and whether they are as strong 
as the ones on planar 3-point functions. We hope to return to these interesting 
issues in the future. 
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A Reflection Positivity 

In this appendix we briefly review the argument for the positivity of 2-point func- 
tions in a theory with a non-real action (as is the case with (2)), and we show 
what positivity implies for the 2-point functions we compute in section 4.1. 

Consider a field theory in (^-dimensional Euclidean space, with fields collectively 
denoted by $ and an action S. Suppose we have a "parity" transformation V that 
satisfies V 2 = 1 and splits the space into two components, with a boundary that is 
given by x = Vx. Let us further assume that V leaves the path integral measure 
X?$ invariant, and that it acts as complex conjugation on the action, in the sense 
that 

S(V$) = S*($). (60) 

Note that V is not a symmetry if the action is not real. Let us denote the fields, 
restricted to either space component, by $±, and the fields on the boundary by $o- 
Now, suppose we have a local operator O that obeys VO{x) = 0(Vx). We can 
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then show that the 2-point function (OO') is positive when evaluated at separated 
points related by parity. Indeed, 

(0{x)0 ] {Vx)) = Jv^e~ s ^O{x)0^{Vx) 

= jv% J[D<S> + ]^ o e- s ^O(x) JpQ-U^e-W-WiVx 

= Jw> J [V$ + ] $b=$o e- s ^O{x) > , (61) 

where in the last step we made the change of variables $_ — > V<&- = $+. Similarly, 
for an operator satisfying VO(x) = —0(Vx), the correlator (OO') is negative. 

For the action (2), which is purely imaginary, we can define a parity transfor- 
mation which has the desired properties: 

x 3 -x 3 , ip^rfip, ^^^7 3 , Af, ->■ (-1)W M . (62) 

Under this V we have — > and J± — > —J±. Reflection positivity then 
implies that 

(jV)(x)jW(Vx))>0, (63) 
(j| 1} (x) J m (Vx)) = -(J ( +\x)J^\Vx)) > . (64) 

Using the known conformal structures [18] for these correlators, we find that the 
general 2-point functions are given by 

(J(°)(x)J(°)(0)) = ^, (65) 



(jW(x)jW(0))=r 1 ^ J X " X \ (66) 



with positive tq and T\. It is easy to verify that this conclusion is consistent with 
our momentum-space results (33) and (37). 
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